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METHOD OF IDENTIFYING STRATEGIES FOR TREATMENT OR 
PREVENTION OF VENTRICULAR FIBRILLATION AND VENTRICULAR 

TACHYCARDIA 

[0001] This application is a national stage application under 35 U.S.C. § 371 

from PCT Application No. PCT/US2005/002848, filed February 2, 2005, which 
claims priority benefit of U.S. Provisional Patent Application Serial No. 60/541,531, 
filed February 3, 2004, which is hereby incorporated by reference in its entirety. 
[0002] This application was made, at least in part, with funding received from 

the National Institutes of Health under Contract No. HL62543. The U.S. government 
has certain rights in this invention. 

FIELD OF THE INVENTION 

[0003] The present invention relates to a method of identifying strategies for 

the treatment or prevention of ventricular fibrillation and ventricular tachycardia. 

BACKGROUND OF THE INVENTION 

[0004] Catastrophic heart rhythm disorders are among the leading causes of 

death in the United States. The most dangerous of these arrhythmias is ventricular 
fibrillation, a disturbance in which disordered wave propagation causes a fatal 
disruption of the synchronous contraction of the ventricle. Although the exact 
mechanism for fibrillation is still being debated, one theory proposes that fibrillation 
is a state of spatiotemporal chaos consisting of the perpetual nucleation and 
disintegration of spiral waves (Chen et al., "Mechanism of Ventricular Vulnerability 
to Single Premature Stimuli in Open-Chest Dogs," Circ Res 62:1991-209 (1988) and 
Witkowski et al., "Spatiotemporal Evolution of Ventricular Fibrillation," Nature 
392:78-82 (1998)), in association with a period doubling bifurcation of local electrical 
properties (Garfinkel et al., "Preventing Ventricular Fibrillation by Flattening Cardiac 
Restitution," Proc Natl Acad Sci USA 97:6061-6 (2000); Gilmour et al., "Electrical 
Restitution, Critical Mass, and the Riddle of Fibrillation," J Cardiovasc 
Electrophysiol 10:1087-9 (1999); Karma, A., "Electrical Alternans and Spiral Wave 
Breakup in Cardiac Tissue," Chaos 4:461-72 (1994); and Panfilov, A., "Spiral 
Breakup as a Model of Ventricular Fibration," Chaos 8:57-64 (1998)). Nucleation of 
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the initiating spiral wave pair is caused by local conduction block (wave break) 
secondary to spatial heterogeneity of refractoriness in the ventricle (Arce et al., 
"Alternans and Higher-Order Rhythms in an Ionic Model of a Sheet of Ischemic 
Ventricular Muscle;' Chaos 10:41 1-26 (2000); Chen et aL, "Mechanism of 
Ventricular Vulnerability to Single Premature Stimuli in Open-Chest Dogs," Circ Res 
62:1991-209 (1988); Sampson et aL, "Simulation and Prediction of Functional Block 
in the Presence of Structural and Ionic Heterogeneity," Am J Physiol 281 :H2597-603 
(2001); Winfree, When Time Breaks Down (Princeton, NJ: Princeton University 
Press) (1987); and Witkowski et al., "Spatiotemporal Evolution of Ventricular 
Fibrillation," Nature 392:78-82 (1998)). Until recently, spatial heterogeneity was 
thought to result solely from regional variations of intrinsic cellular electrical 
properties (Sampson et aL, "Simulation and Prediction of Functional Block in the 
Presence of Structural and Ionic Heterogeneity," Am J Physiol 281:H2597-603 

(2001) ; and Yan et aL, "Characteristics and Distribution of M Cells in Arterially 
Perfused Canine Left Ventricular Wedge Preparations," Circulation 98:1921-7 
(1998)) or from stimulation at more than one spatial location (Qu et aL, "Mechanisms 
of Discordant Alternans and Induction of Reentry in Simulated Cardiac Tissue," 
Circulation 102:1664-70 (2000); Watanabe et aL, "Mechanisms for Discordant 
Alternans," J Cardiovasc Electrophysiol 12:196-206 (2001); and Winfree, When Time 
Breaks Down (Princeton, NJ: Princeton University Press) (1987)). However, it is now 
appreciated that purely dynamical heterogeneity can be sufficient to cause conduction 
block during single-site stimulation in both homogeneous one-dimensional models of 
canine heart tissue and in rapidly paced canine Purkinje fibres (Echebarria et aL, 
"Instability and Spatiotemporal Dynamics of Alternans in Paced Cardiac Tissue," 
Phys Rev Lett 88:208101 (2002) and Fox et al, "Spatiotemporal Transition to 
Conduction Block in Canine Ventricle," Circ Res. 90:289-96 (2002)). A similar 
mechanism has been shown to precipitate conduction block and spiral break-up in 
models of homogeneous two-dimensional tissue (Fenton et aL, "Multiple Mechanisms 
of Spiral Wave Breakup in a Model of Cardiac Electrical Activity," Chaos 12:852-92 

(2002) ). 

[0005] The period doubling bifurcation implicated in the transition to 

conduction block is manifest as alternans, a beat-to-beat long-short alternation in the 
duration of the cardiac action potential (Garfinkel et aL, "Preventing Ventricular 
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Fibrillation by Flattening Cardiac Restitution," Proc Natl Acad Sci USA 97:6061-6 
(2000); Gilmour et aL, "Electrical Restitution, Critical Mass, and the Riddle of 
Fibrillation/' J Cardiovasc Electrophysiol 10:1087-9 (1999); Karma, A., "Electrical 
Alternans and Spiral Wave Breakup in Cardiac Tissue," Chaos 4:461-72 (1994); 
Panfilov, A., "Spiral Breakup as a Model of Ventricular Fibration," Chaos 8:57-64 
(1998); and Watanabe et aL, "Mechanisms for Discordant Alternans," J Cardiovasc 
Electrophysiol 12:196-206 (2001)). Previous investigators have hypothesized that 
alternans can be accounted for by a simple uni-dimensional return map called the 
action potential duration restitution function (Chialvo et al., "Low Dimensional Chaos 
in Cardiac Tissue," Nature 343:653-7 (1990); Chialvo et aL, "Non-Linear Dynamics 
of Cardiac Excitation and Impulse Propagation," Nature 330:749-52 (1987); Guevara 
et aL, ''Electrical Alternans and Period Doubling Bifurcations," IEEE Comput Cardiol 
562:167-70 (1984); and Nolasco et aL, "A Graphic Method for the Study of 
Alternation in Cardiac Action Potentials,' ' J Appl Physiol 25:191-6 (1968)). This 
hypothesis assumes the duration D of an action potential depends only on its 
preceding rest interval / through some function f (I) that is measured experimentally. 
If the D restitution function has a slope >1 , then a period doubling bifurcation occurs 
for some value of the stimulus period 1\ where T= D + I. The velocity Fat which an 
action potential propagates can also be described by a restitution function, where V = 
c(I). 

[0006] It has also been shown previously that the combination of a steeply 

sloped action potential duration ("APD") restitution function and a monotonically 
increasing conduction velocity ("CV") restitution function is sufficient to produce 
dynamical conduction block during sustained pacing at a short cycle length (Fox et al, 
''Conduction Block in One-dimensional heart Fibers," Phys. Rev. Lett. 89:198101 
(2002); Fox et al, "Spatiotemporal Transition to Conduction Block in Canine 
Ventricle," Circ Res, 90:289-96 (2002)). This observation may provide a generic 
mechanism for wave break and the onset of ventricular tachycardia and fibrillation. 
However, it is unlikely that the conditions used to demonstrate this phenomenon 
experimentally apply to the clinical situation, where the induction of ventricular 
tachyarrhythmias typically is associated with the interruption of normal cardiac 
rhythm by only a few premature beats. A single premature beat is sufficient to cause 
spatial heterogeneity in the form of discordant alternans (Watanabe et aL, 
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"Mechanisms for Discordant Alternans," J Cardiovasc Electrophysiol 12:196-206 
(2001)), but the conditions required for the development of conduction block in this 
setting have not been studied extensively. 

[0007] Other studies, building on earlier theoretical work by Krinsky, Winfree 

and colleagues (Krinsky et al, "Votices with Linear Cores in Mathematical Models of 
Excitable Media," Physica A 188:55-60 (1992) and Winfree, A., "Evolving 
Perspectives During 12 Years of Electrical Turbulence," Chaos 8:1-20 (1998)) and 
experiments by Allessie (Allessie et al., "Circus Movement in Rabbit Atrial Muscle as 
a Mechanism of Tachycardia. III. The 'Leading Circle' Concept: a New Model of 
Circus Movement in Cardiac Tissue Without the Involvement of an Anatomical 
Obstacle," Circ Res 41 :9-18 (1977)), have suggested that spiral wave re-entry could 
be the 'engine' that drives ventricular fibrillation ("VF") (Frazier et al., "Stimulus- 
Induced Critical Point. Mechanism for Electrical Initiation of Reentry in Normal 
Canine Myocardium," J Clin Invest 83:1039-1052 (1989); Witkowski et ah, 
"Spatiotemporal Evolution of Ventricular Fibrillation," Nature 392:78-82 (1998); 
Weiss et al, "Chaos and the Transition to Ventricular Fibrillation: a New Approach 
to Antiarrhythmic Drug Evaluation," Circulation 99:2819-2826 (1999); Chen et al., 
"Mechanism of Ventricular Vulnerability to Single Premature Stimuli in Open-Chest 
Dogs," Circ Res 62:1 191-1209 (1988); Gilmour et al., "Electrical Restitution, Critical 
Mass, and the Riddle of Fibrillation," J Cardiovasc Electrophysiol 1 0: 1087-1 089 
(1999); and Pertsov et al., "Spiral Waves of Excitation Underlie Reentrant Activity in 
Isolated Cardiac Muscle," Circ Res 72:631-650 (1993)). Although there is substantial 
evidence that spiral wave re-entry contributes significantly to the induction and 
maintenance of VF, the exact mechanisms by which spiral waves sustain VF is 
currently being debated. 

[0008] The present invention is directed towards correctly identifying the 

mechanisms (or more likely, mechanisms) by which spiral waves cause VF and the 
development of pharmacological approaches to VF treatment and prevention. 



SUMMARY OF THE INVENTION 

[0009] The present invention relates to a method of evaluating the effect of a 

physiological condition on the occurrence of ventricular fibrillation. This method 
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involves providing a test system and initiating a ventricular fibrillation inducing 
sequence in the test system by interrupting normal sinus heart rhythm with premature 
electrical stimuli. Following the initiation of ventricular fibrillation a ventricular 
recovery sequence is initiated. The test system is subjected to a physiological 
condition before, during, and after initiating the ventricular fibrillation recovery 
sequence and the physiological conditions which affect ventricular fibrillation 
recovery are identified. 

[0010] Another aspect of the present invention relates to method of identifying 

treatment candidates as therapeutic strategies for the prevention of ventricular 
fibrillation. This method involves providing a test system and initiating a ventricular 
fibrillation inducing sequence in the test system by interrupting normal sinus heart 
rhythm with premature electrical stimuli before, during, and after administering the 
treatment candidate to the test system. Treatment candidates which prevent the 
initiation step from inducing ventricular fibrillation are identified as therapeutic 
strategies for prevention of ventricular fibrillation. 

[0011] A further aspect of the present invention relates to a method of 

identifying treatment candidates as therapeutic strategies for treating ventricular 
fibrillation. This method includes providing a test system and initiating a ventricular 
fibrillation inducing sequence in the test system by interrupting normal sinus heart 
rhythm with premature electrical stimuli, resulting in ventricular fibrillation in the test 
system. Treatment candidates are administered to the test system undergoing 
ventricular fibrillation and treatment candidates which modulate ventricular 
fibrillation are identified as therapeutic strategies for treatment of ventricular 
fibrillation. 

[0012] The present invention also relates to a method of evaluating the 

predisposition of a subject for the induction of ventricular fibrillation from a condition 
of ventricular tachycardia. This method involves providing a subject in ventricular 
tachycardia and monitoring the electrical stimuli in the heart of the subject. It is then 
determined if a sequence of rest interval values correspond to rest interval values 
predicted to lead to ventricular fibrillation consistent with the histogram of Figure 7. 
[0013] Another aspect of the present invention relates to a method of 

identifying treatment candidates as therapeutic strategies for preventing ventricular 
tachycardia from developing into ventricular fibrillation. This involves providing a 
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subject in ventricular tachycardia and monitoring the electrical stimuli in the heart of 
the subject. It is then determined if the initial 3 stimuli in groups of 4 stimuli 
correspond to rest interval values predicted to lead to ventricular fibrillation consistent 
with the histogram of Figure 7. Treatment candidates which prevent occurrence of a 
fourth stimuli corresponding to a rest interval value predicted to lead to ventricular 
fibrillation consistent with the histogram of Figure 7 are identified. Treatment 
candidates which prevent ventricular tachycardia from becoming ventricular 
fibrillation are identified as therapeutic strategies for preventing ventricular 
fibrillation. 

[0014] One of the hypotheses to account for the apparent link between spiral 

waves and VF is the restitution hypothesis, which proposes that VF is caused by the 
breakup of a single spiral wave into multiple self-perpetuating wavelets (Weiss et al., 
"Chaos and the Transition to Ventricular Fibrillation: a New Approach to 
Antiarrhythmic Drug Evaluation," Circulation 99:2819-2826 (1999); Gilmour et al., 
"Electrical Restitution, Critical Mass, and the Riddle of Fibrillation/' J Cardiovasc 
Electrophysiol 10:1087-1089 (1999); Garfinkel et al., "Preventing Ventricular 
Fibrillation by Flattening Cardiac Restitution," Proc Natl Acad Sci USA 97:6061- 
6066 (2000); Riccio et al., "Electrical Restitution and Spatiotemporal Organization 
During Ventricular Fibrillation," Circ Res 84:955-963 (1999); Panfilov et al., 
"Ventricular Fibrillation: Evolution of the Multiple Wave Hypothesis," Philos Trans 
R SocLondSerA 359:1315-1325 (2001); Fenton et al., "Multiple Mechanisms of 
Spiral Wave Breakup in a Model of Cardiac Electrical Activity," Chaos 12:852-892 
(2002); Karma, A., "Spiral Breakup in Model Equations of Action Potential 
Propagation in Cardiac Tissue," Physc Rev Lett 71:1 103-1 106 (1993); and Karma, A., 
"Electrical Alternans and Spiral Wave Breakup in Cardiac Tissue," Chaos 4:461-472 
(1994), which are hereby incorporated by reference in their entirety). This 
mechanism is similar to that proposed decades ago by Moe (Moe, G., "A Computer 
Model of Atrial Fibrillation," Am Heart J 67:200-220 (1964), which is hereby 
incorporated by reference in its entirety), but with one important difference: it may 
occur in intrinsically homogeneous cardiac tissue. The transition from normal planar 
wave excitation to a single spiral wave and ultimately to multiple wavelets is thought 
to underlie the transition from normal sinus rhythm to ventricular tachycardia and VF 
characteristic of patients who succumb to sudden death (Figure 1). 
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[0015] The exact mechanism for the breakup of spiral waves is unknown. 

However, there is considerable evidence that breakup is closely related to APD 
restitution, which is the relationship between action potential duration (APD) and 
diastolic interval (DI, the time interval between action potentials), where APD is 
determined by the preceding DI (Figure 2). During pacing of cardiac tissue at 
progressively shorter cycle lengths, APD decreases, until at sufficiently short cycle 
lengths a period doubling bifurcation occurs and APD begins to alternate between a 
long duration and a short duration, a phenomenon known as APD alternans (Figure 2) 
(Moe, G., "A Computer Model of Atrial Fibrillation," Am Heart J 67:200-220 (1964) 
(Guevara et aL, "Electrical Alternans and Period Doubling Bifurcations," IEEE Camp 
Cardiol 562:167-170 (1984); Chialvo et aL, "Low Dimensional Chaos in Cardiac 
Tissue," Nature 343:653-657 (1990); Karagueuzian et aL, "Action Potential Alternans 
and Irregular Dynamics in Quinidine-Intoxicated Ventricular Muscle Cells. 
Implications for Ventricular Proarrhythmia," Circulation 87:1661-1672 (1993); 
Nolasco et aL, "A Graphic Method for the Study of Alternation in Cardiac Action 
Potentials," JAppl Physiol 25:191-196 (1968); Rosenbaum et aL, "Electrical 
Alternans and Vulnerability to Ventricular Arrhythmias," New Engl J Med 330:235- 
241 (1994); Pastore et aL, "Mechanism Linking T-Wave Alternans to the Genesis of 
Cardiac Fibrillation," Cir Res 99:1385-1394 (1999), which are hereby incorporated by 
reference in their entirety). The development of APD alternans requires that the slope 
of the restitution relation exceed 1, whereas if the slope of the restitution is less than 
1 , alternans will not occur. 

[0016] The slope of the restitution relation, and the corresponding presence or 

absence of APD alternans, has been linked to spiral wave stability. If the slope of the 
APD restitution relation is < 1, a spiral wave tends to stabilize and produce a periodic 
rhythm, the manifestation of which might be monomorphic ventricular tachycardia. 
If, however, the slope of the APD restitution relation is > 1, a single spiral wave may 
disintegrate into many spiral waves, manifest as VF. These transitions are illustrated 
in Figure 3, using results generated by a computer model. Once initiated, a single 
stable spiral wave can be destabilized by increasing the slope of the restitution 
relation from < 1 to > 1 . Similarly, multiple wavelets can be induced to coalesce into 
a single spiral wave by reducing the slope of the restitution relation from > 1 to < L 
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Thus, a reduction of the slope of the restitution relation from > 1 to < 1 is expected to 
prevent the induction of VF and to convert existing VF into a periodic rhythm. 
[0017] Several recent studies have provided experimental and theoretical 

support for a causal relationship between APD restitution and VF (Garfinkel et al., 
"Preventing Ventricular Fibrillation by Flattening Cardiac Restitution," Proc Natl 
AcadSci USA 97:6061-6066 (2000); Riccio et al., "Electrical Restitution and 
Spatiotemporal Organization During Ventricular Fibrillation," Circ Res 84:955-963 

(1999) ; Koller et al., "Dynamic Restitution of Action Potential Duration During 
Electrical Alternans and Ventricular Fibrillation," Am J Physiol 275:H1635-H1642 
(1998); Koller et al., "Effects of [K+] G on Electrical Restitution and Activation 
Dynamics During Ventricular Fibrillation," Am J Physiol Heart Circ Physiol 
279:H2665-H2672 (2000); Fox et ah, "Spatiotemporal Transition to Conduction 
Block in Canine Ventricle," Circ Res. 90:289-296 (2002); Chen et al., "Spirals, 
Chaos, and New Mechanisms of Wave Propagation," Pacing Clin Electrophysiol 
(PAC) 20:414-421 (1997); and Panfilov, A., "Spiral Breakup as a Model of 
Ventricular Fibrillation," Chaos 8:57-64 (1998), which are hereby incorporated by 
reference in their entirety. If the APD restitution relation contains a region of slope > 
1 , APD alternans and VF are induced by pacing at short cycle lengths. If the slope of 
the APD restitution relation is reduced to less than one, APD alternans is suppressed 
and VF is not induced. Furthermore, if a fibrillating ventricle is exposed to a drug 
that reduces the slope of the restitution to less than 1 , VF is converted to a periodic 
rhythm sustained by a single stable spiral wave (Garfinkel et al., "Preventing 
Ventricular Fibrillation by Flattening Cardiac Restitution," Proc Natl Acad Sci USA 
97:6061-6066 (2000); Riccio et al., "Electrical Restitution and Spatiotemporal 
Organization During Ventricular Fibrillation," Circ Res 84:955-963 (1999); and 
Koller et al., "Effects of [K+] 0 on Electrical Restitution and Activation Dynamics 
During Ventricular Fibrillation," Am J Physiol Heart Circ Physiol 279:H2665-H2672 

(2000) , which are hereby incorporated by reference in their entirety). 
[0018] Unfortunately, interventions that suppress VF experimentally 
(verapamil (Riccio et al., "Electrical Restitution and Spatiotemporal Organization 
During Ventricular Fibrillation," Circ Res 84:955-963 (1999), which is hereby 
incorporated by reference in its entirety), bretylium (Garfinkel et al., "Preventing 
Ventricular Fibrillation by Flattening Cardiac Restitution," Proc Natl Acad Sci USA 



9 



97:6061-6066 (2000), which is hereby incorporated by reference in its entirety) and 
hyperkalemia (Koller et al., "Effects of [K+] 0 on Electrical Restitution and Activation 
Dynamics During Ventricular Fibrillation/' Am J Physiol Heart Circ Physiol 
279:H2665-H2672 (2000), which is hereby incorporated by reference in its entirety) 
have unwanted actions that severely limit their clinical utility. Nevertheless, the 
effects of these interventions on restitution and VF have provided valuable insights 
regarding new, potentially more clinically relevant, drug targets. In particular, 
attempts to understand the mechanism by which calcium channel blockers alter APD 
restitution have led to the realization that increasing selected outward repolarizing 
currents also may flatten restitution, as presented in the examples. 
[0019] The physiological mechanism associated with the initiation of 

ventricular fibrillation are broadly debated. The present invention affords the ability 
to characterize and identify the events initiating ventricular fibrillation. As 
demonstrated by the examples, the present invention also affords novel methods of 
identifying candidates for the treatment and prevention of ventricular fibrillation. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0020] Figure 1 shows the proposed evolution of cardiac wave propagation 

patterns underlying the transition from a normal cardiac rhythm to ventricular 
fibrillation. Transition from normal sinus rhythm (NSR) to ventricular tachycardia 
(VT) and ventricular fibrillation (VF), as recorded on the surface electrocardiogram 
are shown. Corresponding transitions from a planar wave to a single wave to multiple 
spiral waves are shown. 

[0021] Figures 2A-C show the restitution of action potential duration (APD). 

Figure 2A shows the definition of APD, diastolic interval (DI), and cycle length (CL). 
The duration of the n+lth action potential (APD n +i) is a function of the preceding DI 
(D )n ). Figure 2B shows the plot of APD versus preceding DI. Figure 2C shows the 
plot of APD versus CL. Crosses are APDs during the CLs in which no alternans 
occurs. During APD alternans, unfilled circles are APDs during the long action 
potential, and filled circles are APDs during the short action potential. Modified, with 
permission, from Koller et al. (Koller et al., "Dynamic Restitution of Action Potential 
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Duration During Electrical Alternans and Ventricular Fibrillation," Am J Physiol 
275:H1635-H1642 (1998), which is hereby incorporated by reference in its entirety). 
[0022] Figures 3 A-B show the wave patterns in a computer model of cardiac 

activation. Sequential snapshots (A-l) of activation during the development of a 
stable single spiral wave under conditions in which the slope of the restitution relation 
for action potential duration is < 1 (Figure 3 A). The transition form a single spiral 
wave into many spiral waves for a restitution slope >1 (A-O, with some intervening 
steps deleted: Figure 3B). 

[0023] Figure 4A shows a schematic representation of the relationships 

between stimulus period (T), action potential duration (D), diastolic interval (I) and 
memory (M). Figure 4B shows an example of the stimulation protocol. The model 
fibre was paced for ten beats at a constant S' 1 -S' 1 interval, after which a series of 
premature stimuli (S'2-S'5) was delivered. R indicates the refractory period of a 
cardiac cell, during which the delivery of a stimulus does not produce a propagated 
response. 

[0024] Figure 5 shows the spatial dispersion of repolarization induced in an 

ionic model of a homogenous 1 -dimensional cable by multiple premature stimuli. 
Note that short duration action potentials at the site of stimulation increase in duration 
as they propagate alone the fiber, whereas long duration action potentials shorten. 
The resultant discordant alternans pattern culminates in conduction block. 
[0025] Figure 6 shows an example of a type II conduction block in the 

coupled maps model. Each panel shows I as a function of cell number for the last S' 
- 500 ms beat (left panel) and for beats S'2 = 207 ms, S'3 = 146 ms, S'4 - 148 ms 
and S'5 = 106 ms. Note the difference in scale for each beat. Conduction block, 
corresponding to I < I m i n , occurs near cell 100 on the S'5 beat and is denoted by two 
horizontal black lines. 

[0026] Figure 7 shows a histogram of the incidence of type II conduction 

block in the wild type model, a = 0.2, x D = 28 and p = 17.408. The height at each 
point in (S'2, S'3) space corresponds to the number of blocks found at beat S'3, S ? 4 
or S'5 for a given (S'2, S'3) pair. Total counts - 67 205. S'2 min and S'3 min are the 
minimum time intervals at points S'2 and S'3, respectively, that allow conduction 
down the fiber, td is the time interval between activations used to adjust dependence 
of the D recovery function on L (3 is the velocity variable which was varied to adjust 
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the value of V at I = I m j n . a is the memory variable (which varied between zero and 
unity) which was used to determine the influence of memory on D. 
[0027] Figures 8A-C show the role of V restitution in type II block. Figure 

8 A shows V restitution for four cases: p = 17.408 (y-axis - second line from top; wild 
type), P - 7.704 (y-axis - third line from top; V(I min ) = 50% of V max ) ? p - 2.2028 (y- 
axis - fourth line from top; V(I min ) = 25% of V max ), and p = 30.236 (y-axis - first line 
from top; V(I min ) - 90% of V max ). Figure 8B shows a histogram for V(I min ) - 50% of 
V max . Total counts = 8879 (cf 67 205 for wild type). Figure 8C shows a histogram 
for V(I min ) = 90% of V max , Total counts - 86 032. There were no counts for V(I mm ) = 
25%ofV max . 

[0028] Figures 9A-D show the role of D restitution in type II block. Figure 

9 A shows steady-state D restitution for four cases: x D = 28 (y-axis - third line from 
top; wild type; max. slope = 1.09), x D = 32 (y-axis - second line from top; max. slope 
= 0.97), tq — 40 (y-axis - first line from top; max. slope = 0.82) and td^ 20 (y-axis - 
fourth line from top; max, slope == 1 .45). Figure 9B shows a histogram for x D = 32. 
Total counts = 39 688. Figure 9C shows a histogram for x D = 40. Total counts - 21 
1 85. Figure 9D shows a histogram for x D = 20. Total counts = 364 301 . 
[0029] Figures 10A-D show the role of M in type II block. Figure 10A shows 

the steady-state D restitution for four cases: a = 0.2 (y-axis - second line from top; 
wild type; max. slope = 1 .09), a = 0.58 (y-axis - third line from top; max. slope = 
1.62), a = .08 (y-axis - fourth line from top; max. slope = 2.65) and a = 0.1 (y-axis - 
first line from top; max. slope = 1.09). Figure 10B shows a histogram for a = 0.58. 
Total counts - 48 387. Figure 1 0C shows a histogram for a = 0.8. Total counts = 43 
715. Figure 10D shows a histogram for a = 0.1 . Total counts = 87 433. 
[0030] Figure 1 1 shows examples of MAP recordings obtained from the right 

(RV) and left (LV) ventricles of a closed-chest anesthetized beagle during delivery of 
multiple premature stimuli (CLvf) to the RV. CL V f were chosen to produce maximum 
dynamically-induced dispersion of repolarization, as predicted from the results of the 
coupled maps model shown in Figure 7. 

[0031] Figure 12 shows the incidence of conduction block in Beagles. The 

incidence of conduction block followed the delivery of 4 premature stimuli in the 
coupled maps memory model (Figure 7), and following the delivery of CLvf in 
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anesthetized beagles (filled circles). CL V Finduced VF in all 8 dogs. In contrast, 
patterns of cycle lengths that fell outside CL V f did not induce VF (unfilled circles). 
[0032] Figure 13 shows the spontaneous non-sustained polymorphic 

ventricular tachycardia in a German Shepherd dog. 

[0033] Figures 1 4 A-C show the effects of verapamil on action potential 

duration (APD) restitution and ventricular fibrillation (VF). Figure 14A shows an 
example of the reduction in the slope of the APD restitution relation after exposure to 
verapamil. Figure 14B shows the corresponding reduction in the magnitude of APD 
alternans (AM) during pacing at different basic cycle lengths (BCL) after exposure to 
verapamil. Figure 14C shows the effects of verapamil on the composite Fast Fourier 
Transform (FFT) of monophasic action potential (MAP) recordings during VF in 
arterially perfused canine ventricle. Verapamil, added after 35 min of VF had elapsed, 
converted the broad FFT spectra during VF to a single frequency. Modified, with 
permission, from Riccio et al. (Riccio et aL, "Electrical Restitution and 
Spatiotemporal Organization During Ventricular Fibrillation," Circ Res 84:955-963 
(1999), which is hereby incorporated by reference in its entirety). 
[0034] Figures 15A-B show the putative mechanism for APD alternans and 

the suppression of alternans by calcium channel blockade. Under control conditions, 
alternation of the length of the diastolic intervals (DI) during APD alternans causes 
alternation of the magnitude of Ica, secondary to incomplete recovery of Ic a from 
inactivation during the short DI (Figure 15 A). Reduction of Ic a decreases APD, 
which prolongs DI, so that DI following each APD is sufficiently long to allow 7 
complete recovery of Ic a (Figure 1 5B). In the absence of alternations of Ic a , APD 
alternans in suppressed. 

|0035] Figures 16A-D show the effects of HERG overexpression on 1^ and 

APD alternans in isolated canine cardiac myocytes. Current traces of Ikf (Figure 16A) 
and action potentials during pacing at CL = 200 msec (Figure 1 6B) were obtained 
from myocytes in the day of isolation (Day 0). Current traces of I Kr (Figure 16C) and 
action potentials during pacing at CL = 200 msec (Figure 16D) were obtained from 
myocytes after 3 days in cell culture and infection with adenovirus (Day 3 + HERG). 
HERG overexpression increased Ij<r and suppressed APD alternans. 
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DETAILED DESCRIPTION OF THE INVENTION 
[0036] The present invention relates to a method of evaluating the effect of a 

physiological condition on the occurrence of ventricular fibrillation. This method 
involves providing a test system and initiating a ventricular fibrillation inducing 
sequence in the test system by interrupting normal sinus heart rhythm with premature 
electrical stimuli. Following the initiation of ventricular fibrillation a ventricular 
recovery sequence is initiated. The test system is subjected to a physiological 
condition before, during, and after initiating the ventricular fibrillation recovery 
sequence and the physiological conditions which affect ventricular fibrillation 
recovery are identified. 

[0037] In carrying out this aspect of the present invention, the physiological 

condition to be tested can be velocity restitution, action potential duration restitution, 
or cardiac memory. 

[0038] The step of initiating a ventricular fibrillation inducing sequence is 

carried out by: (1) initiating 4 premature stimuli under conditions effective to initiate 
an excitatory wave; (2) altering velocity recovery function to slow conduction at 
shorter rest intervals; (3) altering action potential duration recovery function to reduce 
sensitivity of action potential duration to a preceding rest interval; or (4) is carried out 
by increasing cardiac memory. The test system can be a test animal, a tissue, a cell 
culture, or an in vitro system. 

[0039] Another aspect of the present invention relates to a method of 

identifying treatment candidates as therapeutic strategies for the prevention of 
ventricular fibrillation. This method involves providing a test system and initiating a 
ventricular fibrillation inducing sequence in the test system by interrupting normal 
sinus heart rhythm with premature electrical stimuli before, during, and after 
administering the treatment candidate to the test system. Treatment candidates which 
prevent the initiation step from inducing ventricular fibrillation are identified as 
therapeutic strategies for prevention of ventricular fibrillation. 
[0040] In carrying out this aspect of the present invention, the treatment 

candidate is a pharmaceutical compound, such as a calcium channel antagonist. 
Alternatively, the treatment candidate is one or more electrical impulses. 
[0041] The step of identifying identifies treatment candidates which achieve 

velocity restitution values consistent with the histograms of Figures 8B and 8C as 
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therapeutic strategies for prevention of ventricular fibrillation. Alternatively, the step 
of identifying identifies treatment candidates which achieve potential duration 
restitution values consistent with the histograms of Figures 8B and 8C as therapeutic 
strategies for prevention of ventricular fibrillation. Finally, the step of identifying 
identifies treatment candidates which achieve cardiac memory values consistent with 
the histograms of Figures 9B and 9C as therapeutic strategies for prevention of 
ventricular fibrillation, 

[0042] Again, the test system can be a subject (i.e., a test animal or patient), a 

tissue, a cell culture, an in vitro system. The step of initiating a ventricular fibrillation 
inducing sequence can be carried out by initiating 4 premature stimuli under 
conditions effective to initiate an excitatory wave. 

[0043] A further aspect of the present invention relates to a method of 

identifying treatment candidates as therapeutic strategies for treating ventricular 
fibrillation. This method includes providing a test system and initiating a ventricular 
fibrillation inducing sequence in the test system by interrupting normal sinus heart 
rhythm with premature electrical stimuli, resulting in ventricular fibrillation in the test 
system. Treatment candidates are administered to the test system undergoing 
ventricular fibrillation and treatment candidates which modulate ventricular 
fibrillation are identified as therapeutic strategies for treatment of ventricular 
fibrillation. 

[0044] In carrying out this aspect of the present invention, the treatment 

candidate is a pharmaceutical compound, such as a calcium channel antagonist. 
Alternatively, the treatment candidate can be one or more electrical impulses. 
|0045] The step of identifying identifies treatment candidates which achieve 

velocity restitution values consistent with the histograms of Figures 8B and 8C as 
therapeutic strategies for prevention of ventricular fibrillation. Alternatively, the step 
of identifying identifies treatment candidates which achieve potential duration 
restitution values consistent with the histograms of Figures 9B and 9C as therapeutic 
strategies for prevention of ventricular fibrillation. As yet another alternative, the step 
of identifying identifies treatment candidates which achieve cardiac memory values 
consistent with the histograms of Figures 10B and 10C as therapeutic strategies for 
prevention of ventricular fibrillation. 
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[0046] In carrying out this aspect of the present invention, the test system can 

be a subject (i.e. a test animal or patient), a tissue, a cell culture, or an in vitro system. 
The step of initiating a ventricular fibrillation inducing sequence is carried out by 
initiating 4 premature stimuli under conditions effective to initiate an excitatory wave. 
[0047] The present invention also relates to a method of evaluating the 

predisposition of a subject for the induction of ventricular fibrillation from a condition 
of ventricular tachycardia. This method involves providing a subject in ventricular 
tachycardia and monitoring the electrical stimuli in the heart of the test animal. It is 
then determined if a sequence of rest interval values correspond to rest interval values 
predicted to lead to ventricular fibrillation consistent with the histogram of Figure 7. 
[0048] Another aspect of the present invention relates to a method of 

identifying treatment candidates as therapeutic strategies for preventing ventricular 
tachycardia from developing into ventricular fibrillation. This involves providing a 
subject in ventricular tachycardia and monitoring the electrical stimuli in the heart of 
the test animal. It is then determined if the initial 3 stimuli in groups of 4 stimuli 
correspond to rest interval values predicted to lead to ventricular fibrillation consistent 
with the histogram of Figure 7. Treatment candidates which prevent occurrence of a 
fourth stimuli corresponding to a rest interval value predicted to lead to ventricular 
fibrillation consistent with the histogram of Figure 7 are identified. Treatment 
candidates which prevent ventricular tachycardia from becoming ventricular 
fibrillation are identified as therapeutic strategies for preventing ventricular 
fibrillation. 

[0049] In carrying out this aspect of the present invention, the treatment 

candidate is a pharmaceutical compound, such as a calcium channel antagonist. 
Alternatively, the treatment candidate can be one or more electrical impulses. 
[0050] In one embodiment of this aspect of the present invention, the 

identifying step identifies treatment candidates which achieve velocity restitution 
values consistent with the histograms of Figures 8B and 8C as therapeutic strategies 
for prevention of ventricular fibrillation. Alternatively, the step of identifying 
identifies treatment candidates which achieve potential duration restitution values 
consistent with the histograms of Figures 9B and 9C as therapeutic strategies for 
prevention of ventricular fibrillation. A further alternative is that the identifying step 
identifies treatment candidates which achieve cardiac memory values consistent with 
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the histograms of Figures 10B and IOC as therapeutic strategies for prevention of 
ventricular fibrillation. 

[0051] Previous theoretical and experimental studies of the impact of dynamic 

heterogeneity of repolarization on the development of conduction block were 
conducted in computer models of a homogenous 1 -dimensional cable and in isolated 
Purkinje fibers during fixed pacing. As the pacing cycle length was shortened 
progressively, a transition from 1:1 phase locking to a concordant alternans (2:2 
locking) to discordant alternans to intermittent 2:1 conduction block occurred (Fox et 
al, "Spatiotemporal Transition to Conduction Block in Canine Ventricle," Circ Res. 
90:289-96 (2002), which is hereby incorporated by reference in its entirety). These 
results confirmed the predictions of previous computer modeling studies of these 
phenomena (Watanabe et ak, "Mechanisms for Discordant Alternans," J Cardiovasc 
Electrophysiol 12:196-206 (2001), which is hereby incorporated by reference in its 
entirety) and extended experimental observations that the development of discordant 
alternans was associated with local conduction block and the initiation of reentry 
(Pastore et ak, "Mechanism Linking T-Wave Alternans to the Genesis of Cardiac 
Fibrillation," Circulation 99:1385-94 (1999); Pastore et ak, "Role of Structural 
Barriers in the Mechanism of Alternans-Induced Reentry," Circ Res. 87:1 157-63 
(2000), which are hereby incorporated by reference in their entirety). In addition, 
they demonstrated that complex dynamics was possible in a homogenous 1- 
dimensional cable, in the absence of intrinsic heterogeneity or anisotropy, and that the 
mechanism for the dynamics could be described using simple unidimensional maps 
for the restitution of APD and CV. 

[0052] The evolution of dynamics during fixed pacing occurs relatively 

gradually, which allows this phenomenon to be studied in the quasi-steady-state, As a 
result, a detailed description of the dynamics and identification of the underlying 
mechanism is possible. It is unlikely, however, that prolonged activation at a short 
constant cycle length initiates VF in patients at risk for sudden death. Therefore, 
investigations into a scenario that is more common in such patients - initiation of VF 
following a run of non-sustained VT, were begun. In the present invention, computer 
models of a homogenous 1 -dimensional cable were generated using a coupled maps 
memory model (Fox et ak, "Period-Doubling Instability and Memory in Cardiac 
Tissue," Phys Rev Lett 89:138101 (2002), which is hereby incorporated by reference 
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in its entirety) and an ionic model (Fox et al., "Ionic Mechanism of Electrical 
Alternans," Am J Physiol Heart Circ Physiol 282 :H5 16-30 (2002), which is hereby 
incorporated by reference in its entirety). The models were paced at a cycle length of 
1000 or 500 ms for 50 beats, after which a series of premature stimuli was delivered 
(Figures 4 and 5). Premature stimuli delivered at short DI (e.g., S 9 2 in Figure 5) 
elicited action potentials at the site of stimulation that, in accordance with the relevant 
APD and CV restitution relations were short in duration and conducted slowly. Slow 
conduction along the fiber provided progressively more recovery time as the action 
potential conducted down the fiber. Accordingly, DI increased along the fiber, as did 
the corresponding APD, resulting in a discordance between APD at the site of 
stimulation and at the distal end of the fiber. Conversely, premature stimuli delivered 
at longer DI (S^ in Figure 5) produced longer duration action potentials that 
conducted more rapidly. The more rapidly conducting action potentials impinged 
progressively on the preceding DI as they conducted along the fiber, creating a shorter 
DI and corresponding APD. Once again, discordant APD between the site of 
stimulation and the distal end of the fiber resulted. The spatial dispersion of 
repolarization caused by this sequence of cycle lengths could be amplified by the 
delivery of multiple premature stimuli, eventually resulting in conduction block (S'5 
in Figure 5). 

[0053] To more completely characterize this sequence of events, all possible 

combinations of SVS's that produced conduction block in the coupled maps model 
(block also occurred with many fewer combinations of S VS'3 and many more 
combinations of SVS'7) were determined. Although approximately 67,000 
combinations produced conduction block, the distribution of such combinations was 
not random, but was clustered according to specific S'2 and S'3 intervals, as shown in 
Figure 7. The incidence of block was greatest when S'2 was delivered at a coupling 
interval near S'lS^mm (the shortest S'iS^ interval that produced a propagated 
response) and S'3 was delivered at a coupling interval approximately 50 ms longer 
than S^S^min- Although not shown on this 3-dimensional plot, the S^SU interval 
most often associated with block was near S^S^min, whereas the S^S's interval was 
approximately 10 ms longer than S^S'smin. This short-long-short-intermediate pattern 
of stimulation produced discordant alternans and conduction produced discordant 
alternans and conduction block of the type shown in Figure 5. Also determined was 
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the impact of changes in the slope of the APD restitution relation, the magnitude of 
memory, the shape of the CV restitution relation and the strength of cell coupling. In 
general, the incidence of conduction block was reduced by decreasing the slope of the 
APD restitution relation, increasing memory, increasing cell coupling and extending 
the CV restitution relation to permit slow conduction (thereby preventing "head-tail" 
interactions) (Fox et al., Dynamic Mechanism for Conduction Block in Heart Tissue/' 
New J Phys. 5:101.1-101.14 (2003), which is hereby incorporated by reference in its 
entirety). 

EXAMPLES 

[0054] The following examples are provided to illustrate embodiments of the 

present invention but are by no means intended to limit its scope. 

Example 1 - Methods of Examining the Dynamic Mechanism for Conduction 
Block in Heart Tissue 

[0055] The mechanism for conduction block following the delivery of 

multiple premature stimuli using a one-dimensional computer model of a canine heart 
fiber is presented. The fiber was stimulated at one end at a cycle length near normal 
canine sinus rhythm (S*l = 500 ms). After ten beats at S'l , four premature stimuli 
were delivered (S'2, S'3, S'4 and S'5). The coupling intervals between these 
premature intervals were varied and conduction of the resultant action potentials was 
observed (Figure 4). For each combination of premature stimuli, one of three possible 
outcomes occurred: 1) the stimulus elicited an action potential that propagated down 
the entire fiber; 2) the stimulus did not produce an action potential (Type I block); 3) 
the stimulus elicited an action potential that blocked before reaching the end of the 
fiber (Type II block). 

[0056] The first two cases would not be conducive to the development of 

wave break and initiation of reentry, in that conduction of the premature response 
either does not occur at all or occurs equally well everywhere along the fiber. The 
third case, however, could lead to wave break and spiral wave initiation if it occurred 
in two- and three-dimensional tissue, provided the block was local. The development 
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of local block would be facilitated in intact myocardium by twist anisotropy and 
intrinsic heterogeneity. 

[0057] To assess the vulnerability of the simulated tissue to Type II block, the 

S'2-S 5 3-S'4-S'5 combinations were applied in the following way: the S' l-S'2 
interval was varied from the minimum value that conducted (S'l-S'2 min ) to S' 1- 
S'2 m in+ 20 ms. For each S'l-S'2 interval, the S'2-S^3 interval and the S'3-S'4 
interval were varied in combination from the minimum value that conducted up to a 
value of 250 ms for each interval. Finally, for each S'2-S'3-S'4 combination that 
conducted, the S'4-S'5 interval was varied from the minimum that generated an 
action potential at the site of stimulation to the minimum value that conducted down 
the entire fiber. If no S'2-S'3-S ? 4-S'5 combination was found that produced 
conduction block for more than a 20 ms window in either the S'2-S'3 interval or the 
S'3-S'4 interval, the search was halted for that interval. All intervals were varied in 
steps of 1 ms. 

Example 2 - Coupled Maps Model 

[0058] One aspect of the present invention is conducted using a coupled maps 

model of a one-dimensional cardiac fiber that has been described in detail elsewhere ( 
Fox et aL, "Period-Doubling Instability and Memory in Cardiac Tissue," Phys Rev 
Lett 89:138101 (2002) and Fox et al, "Spatiotemporal Transition to Conduction Block 
in Canine Ventricle," Circ Res, 90:289-96 (2002), which are hereby incorporated by 
reference in their entirety). Briefly, the model is based on the equation 

I^ix^T^x^-D^ix,). (1) 



T n +i (Xi) was the time interval between activations of site x h It was determined by 
including the time delays caused by the propagation from the pacing site to site , 
which yielded 
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x was the time interval between activations applied to the pacing site and Ax^O.l was 
the length of a single cell (time units in ms and space units in mm). The conduction 
velocity V n (x$ depended only on / through the velocity recovery function V„ =■■ cflj 
given by c(I) = V mwc (l-exp(-(I+fi)/S)). V max ^0, 72, 5=14, and jB was varied to adjust the 
value of V at I=I m m. D was determined locally based on a memory model mapping (8) 
given by: 



T m was the time constant of accumulation and dissipation of memory (both constants 
were chosen to be the same). A = 88, B = 122, C = 40, and r m = 180 (Fox et al. 3 
"Period-Doubling Instability and Memory in Cardiac Tissue," Phys Rev Lett 
89:138101 (2002), which is hereby incorporated by reference in its entirety), a and 
td were varied in this study, a (which varied between 0 and 1) determined the 
influence of memory on D, and x D was used to adjust the dependence of the D 
recovery function on / (Fox et al, "Period-Doubling Instability and Memory in 
Cardiac Tissue," Phys Rev Lett 89:138101 (2002), which is hereby incorporated by 
reference in its entirety). 

[0059] Coupling between sites was included by using the diffusion terms from 

Echebarria and Karma (Echebarria et al., "Instability and Spatiotemporal Dynamics of 
Alternans in Paced Cardiac Tissue," Phys Rev Lett 88:208101 (2002), which is hereby 
incorporated by reference in its entirety). These terms modeled the electrotonic 
current that flowed out of (into) a cell and into (out of) its neighbor if the action 
potential of the first cell was longer (shorter) than that of its neighbor. Including 
diffusion then yielded 




(3) 



=f(M„ +i J n ) + t 2 V 2 D n+] -wVD n+1 , 



(4) 
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with ^=1.0 and Discretizing the derivatives produced a tri-diagonal linear 

system of equations that could then be solved easily. The defining equation for the 
model was therefore 

M, 7+1 =g(M„ 9 I n9 D n ) 

D„ +x =f(M n+] J n ) + <* 2 V*D n+] -wVD n+] (5) 
iA Ax A Ax 

7^c<7„ +l (x y )) 7=oc(/„O y )) 

Finally, conduction block was modeled by setting/ = 0 for / < l min = 2. 

Example 3 - Characterization of Type II Conduction Block in Wild-type Model 

[0060] Figure 6 shows a typical example of type II conduction block in the 

baseline model (hereafter called the wild type model). Despite the homogeneous state 
that existed at the end of the S'l stimuli, dynamical heterogeneity developed 
following delivery of S'2, secondary to V restitution. Heterogeneity was manifest as a 
short-to-long gradient in /, which was magnified by subsequent premature stimuli 
until the action potential generated by the S'5 stimulus encountered a region where / 
< Imm and conduction failed. 

[0061] Figure 7 shows a histogram of type II conduction block at various S'2 

and S'3 intervals generated from the wild type model. Note the presence of a large 
peak centered at an S'2 very close to the minimum value for conduction and at an S'3 
approximately 50 ms longer than the minimum value for conduction. Surrounding the 
peak is a flat plateau region of five to ten blocks per bin that extends throughout the 
entire (S'2, S'3) region that was explored numerically. An overwhelming majority of 
these instances of block occurred after the S'5 stimulus. However, six examples of 
'early block' were found after an S'4 stimulus, and one was found after an S'3 
stimulus. The latter explain the presence of the plateau in the histogram; the plateau is 
made up of 'degenerate' counts in which long S'2 or S'3 intervals that maintained the 
initial S"l homogeneity were followed by early block intervals leading to block at 
S'5. 
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[0062] Since conduction block occurs when / < Imm, changes in Imm will 

produce significant changes in the incidence of conduction block. Because Imm is a 
component of the functions that control V restitution, D restitution and M, changes in 
Imm influence all three dynamic variables, Imm was fixed and individually the roles that 
V restitution, D restitution and M played in determining the risk of type II block, as 
assessed by the S'2-S'3-S'4-S'5 protocol, was examined. 

Example 4 - Contribution of V Restitution to Conduction Block 

[0063] V restitution played two important roles in the development of 

conduction block. First, as mentioned above, V restitution was the source of the initial 
heterogeneity after delivery of the S'2 stimulus. The S'2 stimulus was delivered just 
after the refractory period of the last S' 1 beat so that / at the site of stimulation was 
short (see Figure 4), The short / produced slow V, which increased / for cells distal to 
the stimulus site. This initial dispersion in / following S'2 was magnified by steep D 
restitution at each subsequent beat, ultimately leading to conduction block several 
beats later. 

[0064] Second, the value of the V restitution curve at Imm was an important 

parameter in determining the likelihood of conduction block, as illustrated by Figure 
8. Figure 8 A shows V restitution curves with four different values for the V(Imm) (90, 
75, 50, and 25% of Kmax). V(Imm) in the wild type model was 75% of FW Figure 8B 
shows a histogram for V(Imm) = 50% of Fmax, and Figure 8C shows the histogram for 
V(Im\n) = 90%> of fW For V{Imm) = 90% of Kmax, the model had roughly 30% more 
instances of conduction block than in the wild type model, whereas the 
model with the lower cut-off values produced dramatically fewer instances of 
conduction block. In fact, the model with V(Imm) = 25% of Vmax produced no cases of 
type II block. This result did not depend on the steepness of the V restitution curve. 
For example, no cases of type 11 block were found for a model with (5 = -1.4 and 5 = 
2.0, which also had V(Imm) = 25% of PW, but had a much steeper slope. The lack of 
type II block when V(Imm) = 25% of Fmaxcan be understood by noting that if the V 
restitution curve approaches zero at /min, then a propagating wave that approaches a 
region with a very small / can slow down to allow / to increase just in front of it, 
permitting continued conduction. 
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Example 5 - Contribution of D Restitution to Conduction Block 

[0065] As mentioned above, D restitution contributed to conduction block by 

magnifying heterogeneity in /, If the D restitution function had weak dependence on /, 
any initial dispersion of / due to V restitution would be lessened on subsequent 
stimulations. Figure 9 shows the effect of altering D restitution slope by increasing 
and decreasing the parameter x D . Figure 9A shows steady state restitution curves for 
four different values of x D . The wild type model was t d = 28. As shown in Figures 9B 
and 9C, increasing td led to a significantly smaller peak in the histogram, as well as 
fewer total instances of conduction block. Conversely, decreasing x D , which led to a 
steeper steady state restitution slope, had significantly more counts (Figure 9D). 
However, it is important to note that even for the td = 40 model, which had a 
maximum steady state slope of 0.82, instances of type II block were found. 

Example 6 - Contribution of M to Conduction Block 

[0066] Because increasing M in the coupled maps model has been shown to 

eliminate steady state alternans, it was hypothesized that increasing M also would lead 
to fewer instances of type II block. Figure 10 illustrates the effect of changing Mon 
the likelihood of conduction block in the 101,10 model by increasing and decreasing 
the parameter a. Figure 10A shows steady-state restitution for four different values of 
a. The wild type model was a = 0.2. In Figure 10B, a = 0,58, which eliminated 
alternans but maintained a slope £ 1. In fact, the steady state restitution slope was 
1.62, even larger than in the wild type model. In Figure 10C, a = 0.8 and the 
maximum steady state restitution slope was 2.65. Both Figures 10B and 10C show 
smaller peaks and fewer total counts than the wild type model. In Figure 10D, M was 
decreased (a = 0. 1) and the maximum steady state slope is 1 .09. The decreased 
memory model shows more instances of conduction block than the wild type model. 
This result did not depend on the magnitude of D. For example, shifting the curve 
from Figure 10C by adding 70 ms to the function/ (Mn+\, In) did not produce more 
instances of type II block (in fact slightly fewer (39 229) were found). Similarly, 
shifting the curve from Figure 10A down by subtracting 70 ms from/ (Mn+\, In) 
produced roughly the same number of blocks (94 151) as in Figure 10D. Two other 
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observations can be made. First, increasing Mdid not seem to be as effective as 
increasing xo in reducing the likelihood of block. Second, the steady-state D 
restitution function is not predictive of the likelihood of block. 

Example 7 - Contribution of Electronic Interactions to Conduction Block 

[0067] Electrotonic effects have been shown to prevent alternans in some 

models (Cytrynbaum et al., "Stability Conditions for the Traveling Pulse: Modifying 
the Restitution Hypothesis," Chaos 12:788-99 (2002), which is hereby incorporated 
by reference in its entirety). Therefore, changes in electrotonic current may be 
expected to have an effect on the incidence of conduction block. The coupled maps 
model allowed for the complete removal of electrotonic interactions by eliminating 
the diffusion terms in equation (4). As shown in Figure 6, the potential for conduction 
block was present when the previous beat generated a short-to-long pattern in /. This 
pattern led to a short-to-long pattern in the following D. In the wild type model, both 
the first and second derivative diffusion terms tended to decrease the long D and 
increase the short D, reducing the amount of heterogeneity and thereby reducing the 
likelihood of type II block. Accordingly, elimination of the diffusion terms was 
expected to increase the incidence of type II block and this expectation was 
confirmed; total counts of type II block were 88 603. 

[0068] In the present invention, type II conduction block following the 

delivery of multiple premature stimuli was caused by the same sequence of events 
that caused conduction block during sustained rapid pacing in previous studies (Fox et 
al, "Spatiotemporal Transition to Conduction Block in Canine Ventricle," Circ Res. 
90:289-96 (2002), which is hereby incorporated by reference in its entirety). 
Stimulation at a short / produced a gradient of increasing / (and, consequently, of D 
and V) along the fibre. The succeeding stimulus, if delivered at an appropriate 
interval, encountered a progressively decreasing / as it propagated down the fibre (i.e., 
the new wavefront encountered the waveback of the previous excitation). Thus, a 
short- long sequence of D at the site of stimulation was associated with the 
development of a long-short sequence of D at the opposite end of the fibre, similar to 
the pattern established during stable discordant alternans (Pastore et al, "Mechanism 
Linking T-Wave Alternans to the Genesis of Cardiac Fibrillation," Circulation 
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99:1385-94 (1999); Qu et al., "Mechanisms of Discordant Alternans and Induction of 
Reentry in Simulated Cardiac Tissue," Circulation 102:1664-70 (2000); and 
Watanabe et al., "Mechanisms for Discordant Alternans," J Car diovasc 
Electrophysiol 12:196-206 (2001), which are hereby incorporated by reference in 
their entirety). When the gradient in / along the fibre was sufficiently steep in the 
ascending direction, the subsequent action potential encountered an / < /mm and 
conduction block occurred. 

[0069] Previous studies have shown that the location of the block can be 

predicted, given the spatial profile of D for the preceding beat (Sampson et ak, 
"Simulation and Prediction of Functional Block in the Presence of Structural and 
Ionic Heterogeneity," Am J Physiol 281 :H2597-603 (2001), which is hereby 
incorporated by reference in its entirety). If this information is not known, however, 
then developing a quantitative theory for the type II block observed with present 
invention becomes difficult, not only because the model consists of complicated 
nonlinear difference equations in both space and time, but also because the 
phenomenon of interest is a transient. Therefore, the usual tools of nonlinear 
dynamics, such as finding steady-state solutions and analyzing their stability, do not 
apply. Still, inspection of the defining equations produces a qualitative understanding 
of the results presented in this study. The mechanism for conduction block relies on 
three determinants. First, a non-constant V restitution function produces spatial 
heterogeneity in /. Next, this heterogeneity is amplified (or at least it does not decay 
rapidly) due to sensitive dependence of D on /. Finally, the V restitution function does 
not drop sufficiently close to zero at very short / , causing the wave to collide with the 
refractory period of the preceding excitation. 

[0070] Preventing the first aspect of type II block could only be done by 

forcing the ^restitution curve to be constant. Even in that case, type II block could 
still occur if there were intrinsic heterogeneity in the fibre. Therefore, the likelihood 
that type II block could be diminished by two different methods was explored. The 
first method was to alter the V restitution function so that conduction slowed at very 
short /, thus preventing encroachment on the refractory period. It is important to note 
that this method is only viable if the decrease in conduction velocity occurs at short /. 
Decreasing maximum conduction velocity is in fact pro-arrhythmic, because it 
decreases the wavelength (Fenton et al,, "Multiple Mechanisms of Spiral Wave 
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Breakup in a Model of Cardiac Electrical Activity," Chaos 12:852-92 (2002) , which 
is hereby incorporated by reference in its entirety). Although selective slowing of 
conduction at short rest intervals may create a greater dispersion of refractoriness after 
a premature beat, the benefits of this intervention seem to outweigh the risks, at least 
under the conditions of this study. Whether a similar strategy will suppress 
conduction block in 2D or 3D tissue with realistic intrinsic heterogeneity remains to 
be determined. 

[0071] The second method for reducing the likelihood of type II block was to 

reduce the sensitivity of D on /. If the system is uni-dimensional, this is done by 
decreasing the D restitution slope. The lower the slope of the D restitution relation, 
the lower the likelihood of spatial dispersion of repolarization secondary to dispersion 
of /. However, the model used in the present invention did not have a simple 
relationship between the slope of the steady-state D restitution function and the 
likelihood of block. Instead, the important quantity is the sensitivity of D to the 
preceding / , given by 

y a/ a/ 8Ma/' 

(6) 



[0072] The quantity y is the amplification factor relating the dispersion in D to 

the dispersion in the preceding /. If y is small, then any initial dispersion in / will 
decay quickly. If y is large, the dispersion of / can persist or even be amplified. It can 
be shown that at short cycle lengths the first term dominates, giving 
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[0073] M* and /*are the steady state values at a given cycle length. This 

equation is of limited quantitative use since the system never reaches steady state (and 
the cycle length is continuously changing). Still, it suggests that the likelihood of 
block can be reduced by either increasing r D or by increasing a. The simulation results 
support this qualitative argument. This result highlights the fact that the dynamics 
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during premature stimulation cannot be predicted from the steady-state D restitution 
function. In fact, the more relevant experimental measure would be the S' l-S'2 D 
restitution function at very short S'2. The slope of this curve at short / would directly 
measure the dispersion of D after a very short S' l-S'2 interval had produced some 
initial dispersion in /. However, relating the S' l-S'2 curve to any quantitative 
prediction about type II block is not straightforward. 

[0074] Programmed electrical stimulation has been used extensively in the 

past to induce ventricular arrhythmias in patients thought to be at risk for the 
development of ventricular fibrillation (Anderson et al., "Clinical Value of Cardiac 
Electrophysiology Studies, Cardiac Electrophysiology. From Cell to Bedside 2 nd edn, 
ed DP. Zipes and J. Jalife (Philadelphia, PA: Saunders) 1 133-50 (1995), which is 
hereby incorporated by reference in its entirety). However, the stimulation protocol 
used in the clinical studies differs significantly from that used here, in that arrhythmia 
induction was attempted using tightly coupled premature stimuli, which produced 
block similar to the type I block. Once the patient's arrhythmia had been induced 
using multiple premature stimuli, the test was repeated after administration of an 
antiarrhythmic drug. In many cases a drug could be found that would suppress the 
patient's inducible arrhythmia. Yet the results of large scale clinical trials eventually 
revealed that patients sent home on such medications were as or more likely to die 
suddenly than patients who did not receive such treatment (Cast, "Preliminary report: 
Effect of Encainide and Flecainide on Mortality in a Randomized Trial of Arrhythmia 
Suppression After a Myocardial Infarction," in New Engl J Med 321 :406-12 (1989) 
and Waldo et al., "Effect of D-Sotalol on Mortality in Patients with Left Ventricular 
Dysfunction After Recent and Remote Myocardial Infarction. The SWORD 
Investigators. Survival With Oral D-Sotalol," Lancet 348:7-12 (1996), which are 
hereby incorporated by reference in their entirety). Clearly, the induction of 
arrhythmias using this particular stimulation protocol was not predictive of the 
development of ventricular fibrillation. 

[0075] The failure of the standard method of programmed stimulation to 

accurately assess vulnerability to ventricular fibrillation may relate to the fact that this 
method was designed primarily to induce conduction block and re-entry by perturbing 
intrinsic heterogeneity. The potential contribution of dynamic heterogeneity to the 
development of conduction block would not be assessed adequately using such a 



28 



protocol. In addition, interventions that were intended to prevent the induction of 
arrhythmias using the standard protocol, such as slowing of conduction or 
prolongation of refractoriness, might be effective in that regard, secondary to a 
reduction in intrinsic heterogeneity, but might not reduce dynamic heterogeneity. 
Consequently, the beneficial effects of such interventions on arrhythmias induced by 
the standard stimulation protocol might not correspond to a reduction in the incidence 
of ventricular fibrillation, if induction of the latter is influenced importantly by 
dynamic heterogeneity. 

[0076] Dynamically induced heterogeneity and conduction block is one of 

many potential mechanisms for spiral break-up in cardiac tissue (see for example 
(Fenton et al., "Multiple Mechanisms of Spiral Wave Breakup in a Model of Cardiac 
Electrical Activity," Chaos 12:852-92 (2002), which is hereby incorporated by 
reference in its entirety)). Because wave propagation in the heart is extremely 
complex and influenced by many factors, it is unlikely that only one mechanism is 
responsible for induction and maintenance of all cardiac arrhythmias. Nevertheless, 
the results of this study suggest several potential interventions that may reduce the 
probability of arrhythmia induction, including flattening the D restitution slope, 
decreasing Kat short / and increasing M There is experimental support for the idea 
that flattening D restitution suppresses ventricular fibrillation (Riccio et al., 
"Electrical Restitution and spatiotemporal Organization During Ventricular 
Fibrillation," Circ Res 84:955-63 (1999), which is hereby incorporated by reference in 
its entirety), but altering Fand M with the objective of preventing ventricular 
fibrillation has not been attempted. Tests of these options may not be likely in the near 
future; however, given that drugs that increase the recovery of sodium channels at 
depolarized potentials are not currently available, and alteration of Mmay not be 
possible until the ionic basis for this phenomenon has been established. Nevertheless, 
judicious alterations of Z), V, and Mmay hold future promise for new and more 
effective therapies for ventricular fibrillation 

Example 8 - Map Recordings of Beagle Right and Left Ventricles 

[0077] More recently, investigations as to whether dynamically-induced 

spatial dispersion of repolarization promotes the development of conduction block 
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and reentrant excitation in settings where underlying intrinsic heterogeneity is present 
were initiated. To that end, the sequence of premature stimulus intervals predicted by 
the computer model to cause conduction block was delivered to the right ventricles of 
closed chest anesthetized normal beagles and MAP recordings were obtained from the 
right and left ventricles. The dogs initially were paced at a cycle length of 400 ms and 
a single premature stimulus (S^) was delivered at progressively shorter S'iS'2 
intervals until S'2 failed to capture (S'lS^min). The S'iS'2 interval was then set to 
S'iS'2min + 40 ms, and an S' 3 stimulus was delivered at progressively shorter S' 2 S ? 3 
intervals until S'3 m j n was encountered. The procedure was then repeated to find S'4min 
and S'smin. Thereafter, a sequence of premature stimuli consisting of S'iS'2min+ 5 ms, 
S^S^min + 50 ms, S^S^min + 5 ms and S^S'smin + 5-10 ms was delivered. As shown 
in Figure 1 1 , the patterns of APD produced by this sequence of premature cycle 
lengths (CLvf) was similar to that produced in the computer models (c.f Figure5). 
Short duration responses (S) at the RV recording site were associated with longer 
duration responses (L) at the LV site and vice versa and the pattern of discordant 
alternans preceded the development of VF. 

[0078] The delivery of CL V p has produced VF in all 8 dogs tested thus far, 

whereas the delivery of multiple premature stimuli at sequentially shorter cycle 
lengths has failed to induce VF in any of the dogs (Figure 12). Although, for obvious 
reasons, not all 67,000 possible combinations of cycle lengths predicted to cause 
conduction block in the model, nor more than a few of the intervals that are not 
expected to cause block were tested, the highly reproducible induction of VF by CLvf 
lying near the peak of the plot shown in Figure 1 1 suggests that stimulus protocols 
designed to maximize dynamically-induced spatial heterogeneity of repolarization 
may reliably promote the induction of VF. The exact nature of the spatial dispersion 
produced by such a protocol will require more extensive mapping than can be 
accomplished by two MAP recordings, but the observation that CLvf produces 
discordant alternans between at least 2 RV and LV recording sites suggests that the 
pattern of APD generated by the computer models may pertain to the intact heart, 
[0079] As a further test of the hypothesis that dynamically-induced spatial 

heterogeneity is arrhythmogenic, it was determined whether a reduction in the slope 
of the APD restitution relation suppresses the development of spatial heterogeneity 
APD and conduction block during premature stimulation, as shown in the computer 
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model (Fox et al., Dynamic Mechanism for Conduction Block in Heart Tissue," New 
J Phys. 5:101.1-101.14 (2003), which is hereby incorporated by reference in its 
entirety). For these studies, the effects of verapamil, which is known to decrease the 
slope of the dynamic restitution relation (Riccio et al. ? "Electrical Restitution and 
Spatiotemporal Organization During Ventricular Fibrillation," Circ Res 84:955-63 
(1999), which is hereby incorporated by reference in its entirety), in the anesthetized 
closed chest dogs described above (n = 8) were tested. MAP duration (MAPD) was 
measured in RV and LV endocardium during programmed stimulation of the RV 
using 4 appropriately timed premature stimuli (SVS's), following a train of 20 stimuli 
at a constant S'iS'i interval of 400 ms. The stimulation protocol incorporated the 
same short-long-short-intermediate sequence shown in Figure 12, adjusted as 
necessary to account for changes in S'imin-S'smjn produced by verapamil. VF 
induction was attempted during control and 30 minutes after i.v. verapamil (0.1, 0.3 or 
1.0 mg/kg/min). 

[0080] VF induction occurred in all dogs during control and after 0. 1 mg/kg 

verapamil and was associated with discordant alternans of MAPD between RV and 
I.V. In contrast, after the two higher verapamil doses MAPD alternans was reduced 
in magnitude and was concordant and VF could not be induced. These effects were 
not associated with significant changes in the slope of the standard (S' ]S' 2 ) MAPD 
restitution relation, which was 0.41 ±0.17 during control and was 0.47 ± 0,28, 0.56 ± 
0.20 and 0.39 ±0.18 after verapamil. Thus, suppression of spatial APD heterogeneity 
by verapamil was associated with an inability to induce VF using premature stimuli. 
However, the doses of verapamil required to produce this effect had significant 
negative dromotropic, chronotropic and inotropic effects: verapamil increased the PR 
interval (from 84.7 ±3.8 during control to 111.8 ± 12.2, 132.0± 15.0 and 139.5 ± 
18.0 ms after verapamil), decreased heart rate (from 139.0 ± 9.7 to 1 10.0 ± 2.1, 101.6 
± 3.8, and 95.2 ± 8.8 bpm) and decreased mean blood pressure (from 1 16.8 ± 7.3 to 
107.4 ±11.4, 108.0 ±6.3 and 101 ,2 ± 9,4 mmHg). These results suggest that 
suppression of VT-induced spatial electrical heterogeneity may be one means of 
preventing VF, but that verapamil is unlikely to be the drug of choice for producing 
such an effect in patients. 

Example 9 - Ventricular Fibrillation in German Shepherd Dogs 
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[0081] Although the induction of VF by CL V f in normal beagles confirms the 

predictions made by the computer models, this result is somewhat expected, in that 
the data used to construct the models was obtained from normal canine myocardium. 
This model has the additional limitation of requiring external stimuli to induce VF, 
since these dogs do not have spontaneous ectopy. To test the predictions from the 
computer model more rigorously, next it was investigated whether CLvf induced VF 
in German shepherd dogs that display non-sustained polymorphic VT and have an 
inherited predisposition to sudden death (Figure 13). As discussed above, one of the 
interesting and unanswered questions with respect to these animals is why they often 
live for months with multiple episodes of rapid, polymorphic VT before succumbing 
to sudden death. Based on the response of the beagles to CLvf, it seems possible that 
the German shepherds, despite having copious amounts of ectopy, do not die sooner, 
because they do not routinely generate the proper sequence of coupling intervals, 
insofar as promoting dynamically-induced spatial dispersion of repolarization is 
concerned. 

[0082] As an initial test of this idea, 5 affected German shepherd dogs were 

anesthetized and instrumented as described above for the beagles. CLvf were 
determined from measurements of S^min-S'smin and predictions from the beagle-based 
computer model and were delivered to the RV or LV. VF was induced in 3 of the 
dogs, but not in the other two, even after several adjustments of the CL V f* Sequences 
of cycle lengths obtained from Holter recordings during multiple episodes of 
polymorphic VT (CLyr) also were delivered to the dogs. None of these sequences 
induced VF. The sequences during CLvt consisted of progressively shorter cycle 
lengths, in contrast to CL V f, which were short-long-short-short. In addition, even 
during longer runs of non-sustained VT, such as the run shown in Figure 13, a CL vf - 
type pattern did not occur. This observation raises the question of why do the German 
shepherds eventually die suddenly? Unfortunately, despite the hundreds of Holter 
recordings obtained since this colony was established, VF has never been captured. 
To remedy this situation, Reveal Plus© (Medtronic) event recorders were implanted 
in severely affected dogs. Thus far, the development of VF in one dog has been 
captured. The sequence of cycle lengths preceding VF was similar to that used to 
precipitate VF in the programmed stimulation studies, but clearly more episodes of 
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spontaneous VF need to be captured and analyzed before definitive conclusions can 
be drawn. 

Example 10 - Calcium Channel Agonists and Ventricular Fibrillation 

[0083] The suppression of VF by verapamil, in association with a reduction in 

the slope of the APD restitution relation, (Figure 14) suggests that the L-type calcium 
current (Ic a ) plays a key role in restitution and in the development of VF (Riccio et al., 
"Electrical Restitution and Spatiotemporal Organization During Ventricular 
Fibrillation," Circ Res 84:955-963 (1999) and Chudin et al., "Intracellular Ca 2 " 
Dynamics and the Stability of Ventricular Tachycardia," Biophys J 77:2930-2941 
(1999), which are hereby incorporated by reference in their entirety). To better define 
these roles, an ionic model was developed for the cardiac ventricular myocyte (CVM) 
based on the Winslow (Winslow et al., "Mechanisms of Altered Excitation- 
Contraction Coupling in Canine Tachycardia-Induced Heart Failure, IL," Circ Res 
84:571-586 (1999), which is hereby incorporated by reference in its entirety) and 
Luo-Rudy models (Luo et al, "A Dynamic Model of the Cardiac Ventricular Action 
Potential, I: Simulation of Ionic Currents and Concentration Changes," Circ Res 
74:1071-1096 (1994), which is hereby incorporated by reference in its entirety) that 
generates physiologically realistic APD alternans over a wide range of pacing cycle 
lengths (Fox et al., "Ionic Mechanism of Cardiac Alternans," Am J Physiol 282:H5 1 6- 
H530 (2002), which is hereby incorporated by reference in its entirety). Using this 
model, the L-type Ca 2+ current was implicated as an important determinant of APD 
alternans, according to the mechanism shown in Figure 1 5. Upon initiation of pacing 
at a short CL, I Ca is fully recovered prior to the first action potential and activates fully 
during the action potential, which results in a long APD. The long action potential 
subsequently is followed by a short diastolic interval, during which I Ca fails to recover 
completely from Ca 2+ -induced inactivation. Because of decreased availability of Ica, 
the duration of the subsequent action potential is shorter, which results in a longer 
succeeding diastolic interval, more complete recovery of Ic a and a long action 
potential duration. This cycle then repeats, eventually establishing a steady-state 
alternans of Ica and APD. 



33 



Example 11 - Potassium Channel Agonists and Ventricular Fibrillation 

[0084] From the mechanism for APD alternans outlined in Figure 1 5, it was 

anticipated that shortening of APD during rapid pacing will prolong DI, which, in 
turn, might provide adequate time for complete recovery of Ic a - If so, reducing the 
magnitude of APD alternans could be accomplished by increasing outward 
repolarizing currents, rather than by decreasing Ic a . This idea has been supported by 
additional computer modeling studies in which APD alternans was suppressed by 
increasing any one of several outward repolarizing currents (e.g., the rapid (IkO and 
slow (Iks) components of the delayed rectifier and the inward rectifier (IkO) (F° x et 
al., "Ionic Mechanism of Cardiac Alternans," ylm J Physiol 282:H516-H530 (2002), 
which is hereby incorporated by reference in its entirety). Although agonists for these 
currents generally are not available, it might be possible to modify currents such as iKr 
and Iks by increasing phosphatidyl inositol bisphosphate (PIP2) levels (Bian et al., 
"HERG K- Channel Activity is Regulated by Changes in Phosphatidyl Inositol 4,5- 
Bisphosphate," Circ Res 89:1 168-1 176 (2001), which is hereby incorporated by 
reference in its entirety) or by altering the phosphorylation state of the channels 
(Kiehn, J., "Regulation of the Cardiac Repolarizing HERG Potassium Channel by 
Protein Kinase A/' Trends Cardiovasc Med 10:205-209 (2001); Heath et al., "Protein 
Kinase C Enhances the Rapidly Activated Delayed Rectifier Potassium Current, l Kr , 
Through a Reduction in C-Type Inactivation in guinea Pig Ventricular Myocytes/' J 
Physiol 522:391-402 (2000); Marx et ah, "Requirement of a Macromolecular 
Signaling Complex for p Adrenegic Receptor Modulation of the KCNQ1-KCNE1 
Potassium Channel," Science 295:496-499 (2002), which are hereby incorporated by 
reference in their entirety), provided the channels can be phosphorylated without 
concomitant phosphorylation of calcium channels or can be upregulated in the 
presence of a calcium channel antagonist, to offset increased I Ca secondary to 
phosphorylation (Reuter et al., 'The Regulation of Calcium Conductance of Cardiac 
Muscle by Adrenaline," J Physiol 246:49-62 (1977), which is hereby incorporated by 
reference in its entirety). 

[0085] Recently, increasing I Kr as an approach to suppressing APD alternans 

has been examined. I Kr plays an important role in cardiac repolarization, increasing to 
a maximum during phase 3 of the action potential, as the channel recovers from 
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inactivation, and then decreasing, as the electrical driving force decreases and as 
deactivation of the channel increases (Sanguinetti et al., 'Two Components of Cardiac 
Delayed Rectifier K+ Current. Differential Sensitivity to Block by Class III 
Antiarrhythmic Agents," J Gen Physiol 96:195-215 (1990); Rocchetti et al., "Rate 
Dependency of Delayed Rectifier Currents During the Guinea Pig Ventricular Action 
Potential," J Physiol 534:721-732 (2001); Gintant, G., "Characterization and 
Functional Consequences of Delayed Rectifier Current Transient Repolarization," Am 
J Physiol 278:h806-H817 (2000); Gintant, G., "Regional Differences in I K Density in 
Canine Left Ventricle: Role of I Ks in Electrical Heterogeneity," Am J Physiol 
268:H604-H613 (1995); and Clay et al., "A Quantitative Description of the E-4031- 
Sensitive Repolarization Current in Rabbit Ventricular Myocytes," Biophys J 
69:1830-1837 (1995), which are hereby incorporated by reference in their entirety). 
Because iRr contributes minimally to the action potential plateau, increasing I Kr may 
have little or no effect on the Ca 2+ transient. Consequently, by increasing Ijcr it may 
be possible to suppress APD alternans without adversely affecting contractility. 
[0086] At present, however, there are no l Kr agonists available to test this 

hypothesis. To circumvent this problem, isolated ventricular myocytes were infected 
with an adenovirus expressing HERG, the gene that encodes the pore-forming domain 
of Ikh to increase HERG protein expression level and the corresponding I Kr current 
(Nuss et al., "Overexpression of a Human Potassium Channel Suppresses Cardiac 
Hyperexcitability in Rabbit Ventricular Myocytes,' V Clin Invest 103:889-896 (1999); 
and Waldo et al., "Effect of D-Sotalol on Mortality in Patients With Left Ventricular 
Dysfunction After Recent and Remote Myocardial Infarction. The SWORD 
Investigators. Survival with Oral D-Sotalol," Lancet 348:7-12 (1996), which are 
hereby incorporated by reference in their entirety). After verifying that Iicr had 
increased post-infection, the myocytes were paced at rapid rates to determine whether 
increasing I Kr by overexpressing HERG suppressed APD alternans. As shown in 
Figure 16, over expression of HERG markedly increased I Kr , as recorded during 
action potential clamp, and suppressed APD alternans during rapid pacing. Over 
expression of HERG did not, however, significantly reduce I C a, suggesting that this 
approach to suppression of alternans need not be accompanied by a reduction in 
contractility. 
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[0087] The current paradigm regarding 1^ and VF is that blocking I^r is 

expected be anti-arrhythmic, secondary to prolongation of APD and refractoriness 
[Waldo et al., "Effect of D-Sotalol on Mortality in Patients With Left Ventricular 
Dysfunction After Recent and Remote Myocardial Infarction. The SWORD 
Investigators. Survival with Oral D-Sotalol. Lancet 348, 7-12 (1996); Hondeghem et 
al., "Class III Antiarrhythmic Agents Have a Lot of Potential But a Long Way to Go. 
Reduced Effectiveness and Dangeers of Reverse Use Dependence," Circulation 
81 :686-690 (1990), which are hereby incorporated by reference in their entirety]. 
However, this approach has not been successful in preventing VF and, moreover, is 
associated with pro-arrhythmia [Waldo et al., "Effect of D-Sotalol on Mortality in 
Patients With Left Ventricular Dysfunction After Recent and Remote Myocardial 
Infarction. The SWORD Investigators. Survival with Oral D-Sotalol. Lancet 348, 7- 
12 (1996), which is hereby incorporated by referenced in its entirety]. Prolongation 
of APD, as reflected by an increased duration of the QT interval on the ECG, and an 
increase in the incidence of ventricular arrhythmias also are associated with the 
administration of numerous non-cardiac drugs that block I Kr [Mohammad et al., 
"Blockage of the HERG Human Cardiac K+ Channel by the Gastrointestinal 
Prokinetic Agent Cisapride," Am. J. Physiol 723:H2534-H2538 (1997); Roy et al„ 
"HERG, a Primary Human Ventricular Target of the Nonsedating Antihistamine 
Terfenadine," Circulation 94:815-823 (1996); Mitcheson et al., "A Structural Basis 
For Drug-Induced Long QT Syndrome," Proc. Natl. Acad ScL U.S.A. 97:12329- 
12333 (2000); Mitcheson et al., "A Structural Basis For Drug-Induced Long QT 
Syndrome," Proc. Natl. Acad Set U.S.A. 97:12329-12333 (2000), which are hereby 
incorporated by reference in their entirety]. Similarly, inherited loss-of- function 
mutations in iKr are accompanied by a prolongation of the QT interval and by an 
increased risk of lethal ventricular tachyarrhythmias such as torsade de pointes 
[Keating et aL, "Molecular Genetic Insights Into Cardiovascular Disease," Science 
272:681-685 (1996); Curran et al., "A Molecular Basis For Cardiac Arrhythmia: 
HERG Mutations Cause Long QT Syndrome," Cell 80:795-803 (1995); Roden et al., 
"A Plethora of Mechanisms in the HERG-Related Long QT Syndrome. Genetics 
Meets Electrophysiology," Cardiovasc. Res. 44:242-246 (1999); Sanguinetti, et al., 
"A Mechanistic Link Between An Inherited and An Acquired Cardiac Arrhythmia: 
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HERG Encodes the I Kr Potassium Channel," Cell 81 :299~307 (1 995), which are 
hereby incorporated by reference in their entirety]. 

[0088] The cardiac arrhythmias associated with inherited or drug-induced 

abnormalities of l& r are thought to be precipitated primarily by bradycardia-induced 
prolongation of repolarization [Keating et aL, "Molecular Genetic Insights Into 
Cardiovascular Disease," Science 272:681-685 (1996); Curran et aL, "A Molecular 
Basis For Cardiac Arrhythmia: HERG Mutations Cause Long QT Syndrome," Cell 
80:795-803 (1995), which are hereby incorporated by reference in their entirety]. The 
potential mechanisms by which I Kr may facilitate the induction of arrhythmias at slow 
heart rates have been studied extensively [Hondeghem et aL, "Class III 
Antiarrhythmic Agents Have a Lot of Potential But a Long Way to Go. Reduced 
Effectiveness and Dangeers of Reverse Use Dependence," Circulation 81:686-690 
(1990); Mitcheson et aL, "A Structural Basis For Drug-Induced Long QT Syndrome," 
Proc. Natl Acad Set U.S.A. 97:12329-12333 (2000); Roden et aL, "A Plethora of 
Mechanisms in the HERG-Related Long QT Syndrome. Genetics Meets 
Electrophysiology," Cardiovasc. Res. 44:242-246 (1999), which are hereby 
incorporated by reference in their entirety]. However, the contribution of Ikf to 
repolarization during tachycardia, which might also be important for the development 
of cardiac arrhythmias, has not been well characterized. The observation that 
reducing I Kr increases the magnitude of APD alternans may provide an additional 
mechanism to account for the proarrhythmic effects of Ikf blockers, in that increased 
alternans magnitude would be expected to destabilize ventricular tachyarrhythmias, 
leading to the development of VF. Conversely, the observation that increasing I Kr 
reduces the magnitude of APD alternans and the slope of the APD restitution relation 
provides a rationale for the development of a new class of compounds, agonists, 
with the expectation that such compounds may have anti-fibrillatory effects. 
[0089] The latter hypothesis relies, however, on the expectation that Iter can be 

increased sufficiently to reduce the slope of the restitution relation without shortening 
action potential duration to such an extent that contractility is impaired, secondary to a 
shortening of the plateau duration and attenuation of Ic a - In addition, shortening of 
APD may lead to a reduction in the wavelength of reentry circuits (where wavelength 
= refractory period x conduction velocity). Reduction of the wavelength could, in 
turn, precipitate wavebreak by other mechanisms (e.g., so-called "head-tail" 
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interactions, where a wavefront encounters a waveback and fragments) [Hund et al., 
"Dynamics of Action Potential Head-Tail Interaction During Reentry in Cardiac 
Tissue: Ionic Mechanisms," Am. J. Physiol 279:H1869-H1879 (2000), which is 
hereby incorporated by reference in its entirety]. These and other potential drawbacks 
to the use of K-channel agonists to suppress VF remain to be evaluated critically. 
[0090] Until recently, therapy for the prevention of sudden cardiac death has 

been based on the presumption that frequent ventricular ectopy, in particular 
ventricular tachycardia, is a prelude to ventricular fibrillation [Myerburg et al., "Life- 
Threatening Ventricular Arrythmias: the Link Between Epidemiology and 
Pathophysilogy," In Cardiac Electrophysiology: From Cell to Bedside (AZipes, D.P. 
and Jalife, J. eds), pp. 723-731, WB Saunders Co., Philadelphia (1995), which is 
hereby incorporated by reference in its entirety]. Accordingly, drugs that suppress 
inducible or spontaneously occurring ventricular tachycardia are expected to prevent 
sudden death. However, recent large-scale clinical trials have indicated that classes of 
drugs that are effective for the suppression of ventricular tachycardia, e.g., Class I and 
Class III antiarrhythmic drugs, do not prevent sudden death and, in fact, may be 
proarrhythmic [ Waldo et al., "Effect of D-Sotalol on Mortality in Patients With Left 
Ventricular Dysfunction After Recent and Remote Myocardial Infarction. The 
SWORD Investigators. Survival with Oral D-Sotalol. Lancet 348, 7-12 (1996); Echt 
et al., "Mrtality and Morbidity in Patients Receiving Encainide, Flecainide or Placebo: 
the Cardiac Arrhythmia Suppression Trial," New Engl J. Med. 324:781-788 (1991), 
which are hereby incorporated by reference in their entirety]. In contrast, other 
classes of drugs that are not particularly effective for the suppression of most forms of 
ventricular tachycardia, such as beta-adrenergic receptor antagonists [Yusuf et al., 
"Overview of Results of Randomized Clinical Trials in Heart Disease: I. Treatments 
Following Myocardial Infarction," J. Am. Med. Assoc. 260:2088-2095 (1988), which 
is hereby incorporated by reference in its entirety] and calcium channel antagonists 
[Held et al., "Impact of Calcium Channel Blockers on Mortality," In Cardiovascular 
Pharmacology and Therapeutics (Singh, B.N. et al., eds), pp. 525-533, Churchill 
Livingston (1994), which is hereby incorporated by reference in its entirety], may 
reduce mortality from sudden death. 

[0091] If a causal relationship between the slope of the APD restitution 

relation and VF is confirmed, such a result could have significant implications for the 
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pharmacological therapy of sudden cardiac death. Drugs that reduce the slope of the 
restitution relation would be expected to prevent the development of VF, but would 
not be expected to suppress ventricular tachycardia, if ventricular tachycardia is 
caused by some variant of spiral wave reentry. Conversely, drugs that do not reduce 
the slope of the restitution relation would not be expected to prevent VF, although 
they might suppress ventricular tachycardia, perhaps via a mechanism that does not 
involve alteration of restitution kinetics (e.g. slowing of conduction or prolongation of 
refractoriness), 

[0092] Given these expectations, the usefulness of existing agents for the 

prevention of VF could be reevaluated in light of their effects on the restitution 
relation and new drugs targeted against restitution could be developed. Naturally, 
such an effort would need to recognize that factors other than APD restitution 
contribute importantly to the induction and maintenance of VF. Furthermore, drugs 
intended to alter restitution selectively might have additional effects that could offset 
their intended effect. Nevertheless, judicious alteration of the APD restitution slope 
appears to be a promising approach to the treatment of VF and, as such, represents a 
potentially fruitful opportunity for drug development. 

[0093] Electrical restitution is a recent addition to factors that play a key role 

in the development of ventricular tachyarrhythmias. Experimental interventions that 
reduce the slope of the restitution relation have been effective in suppressing VF. 
These results encourage further investigation of altering restitution as a means of 
preventing sudden death, with the realization that clinically useful interventions that 
flatten restitution without having untoward effects have yet to be developed. 
[0094] Although the invention has been described in detail for the purposes of 

illustration, it is understood that such detail is solely for that purpose, and variations 
can be made therein by those skilled in the art without departing from the spirit and 
scope of the invention as defined in the claims that follow. 



